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Abstract. - The fluctuations of the vortex density in a turbulent quantum fluid are deduced 
from local second-sound attenuation measurements. These measurements are performed with a 
micromachined open-cavity resonator inserted across a flow of turbulent He-II near 1.6 K. The 
power spectrum of the measured vortex line density is compatible with a (-5/3) power law. The 
physical interpretation, still open, is discussed. 
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Motivation. — Our understanding of the turbulent 
fluctuations in quantum fluid flows relies on a single direct 
measurement [1]. This lack of local experimental data [2] 
contrasts with the situation for counter-flow turbulence, 
which has been explored by a dozen of fluctuations studies 
( [3-5] ) . It also contrasts with the detailed theoretical pic- 
ture of the Kolmogorov-like cascade which emerged over 
the last decades [8] and with the diversity of numerical 
simulations, such as velocity spectrum studies [9-11]. 

In the isolated experiment mentioned above [1], a minia- 
ture total-head tube was configured to operate as an 
anemometer in a He-II co-flow (or "bulk flow"), at 2.08 
and 1.4 if. In a two-fluid model picture, this probe senses 
a combinaison of the velocities of the superfluid and nor- 
mal components. Its spatial and time resolution (typi- 
cally 2 mm and 1 kHz) enabled to resolve one decade of 
power law scaling below the injection scale. This scaling 
was found to overlap with the scaling in the inertial range 
of classical flows. This strong result suggests that the 
largest scales structures in quantum turbulence undergo a 
Kolmogorov-type cascade. 

Second sound probes have a several decades history as 
an efficient tool to measure quantum vortices line den- 
sity (VLD), in particular in turbulent co-flows (for in- 
stance [12-14] and references within) and counter-flows 
(for instance [3-7] and references within). Unfortunately, 
the size of these sensors and their sidewall positionning 
made impossible space and time resolved measurements of 
flow fluctuations. The aim of this paper is to report such 
a local fluctuations measurement from a micromachined 
miniature second sound resonator. This probe therefore 



completes the inertial range characterization with a fully 
superfluid observable, independent of the earlier velocity 
measurements. 



The flow set-up. — Our flow is a He-II loop confined 
in a nearly cylindrical cryostat and continuously powered 
by a centrifugal pump (see figure 1). Turbulence is probed 
in a $ = 23 mm-diameter, 250mm-long brass pipe, lo- 
cated upstream from the pump . Downstream the pump, 
the fluid returns to the pipe inlet flowing along the outside 
of it. On this return path, a 30mm-long 3mm-cell hon- 
eycomb breaks the spin motion generated by the pump. 
As a test, another 20 mm- long honeycomb has once been 
inserted in the pipe inlet without noticeable changes. A 
Pitot tube is located 130 mm before the pipe outlet. It 
provides a measurement of the mean velocity by mean of 
an in-situ capacitive differential pressure gauge. The use- 
full range of velocity is V = 0.3 — 1.3 m/s. At lower veloci- 
ties flow instabilities are detected and at higher velocities, 
typically at 1.5 m/s, a cavitation threshold is encountered. 
From in-situ measurements with a semiconducting minia- 
ture pressure sensors [16], we estimate a typical 35% veloc- 
ity turbulence ratio in the pipe. With velocity V and pipe 
diameter <£, several Reynolds numbers can be defined us- 
ing different denominators. This multiplicity result from 
the extra degrees of freedom of quantum fluids compared 
to classical fluid. To assess the "instability" of the super- 
fluid, a possible denominator is the quantum of circulation 
hz = h/m — 0.997.10 _7 m 2 /s (h is Planck constant and m 
is the mass of the 4-He atom) which is the only available 
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Fig. 2: Left figure : Single silicon element before the probe 
assembling, with contact leads and heater layer (emitter) or 
thermometer layer (receiver). These layers are deposited over 
a 0.8 mm x 0.8 mm surface. Right figure : Side view of the tip 
of the assembled probe. The heater and thermometer are fac- 
ing each other with a 250 \im- gap between their 15/xra-thick 
and 1.5 mm long support. The helium mainstream flows per- 
pendicularly to the right-side figure. 



Fig. 1: Scheme of the flow loop with probes. 

quantum parameter 1 . The corresponding Reynolds num- 
ber Re K = <$>.V/k falls in the 6.10 4 - 3.10 5 range. 

Fluid temperature is regulated near 1.6 which -in the 
two-fluid model- corresponds to 84% of superfluid [15]. 
Regulation heaters are spatially distributed and located 
between the pump and the main honeycomb, near two 
doped germanium thermometers dedicated to regulation 
and monitoring. Cooling is provided by a pumped he- 
lium bath located above the flow. The excellent effective 
thermal conductivity of He-II makes possible an efficient 
regulation. The helium bath's weight pressurizes the flow 
and makes possible cavitation free operation. 

The second sound probe. — The second-sound 
miniature probe is located 80 mm downstream the pipe 
inlet. It consists in an open cavity through which a small 
fraction of helium flows (see figure 2). Both "mirors" of 
the cavity are lb microns thick silicon beams separated 
by a D ~ 250 \im gap. The length and width of both 
beams is 1.5 mm x 1mm and they are facing each other 
with a lateral positioning within a few tenths of mm typ- 
ically. A granular Al film is deposited over a h ~ 0.8 mm 
square area at the tip of one beam. This film is used as 
transition edge superconducting thermometer. Its criti- 
cal temperature is roughly 1.6 K. On the center of the 
transition, its thermal sensitivity is T/R.dR/dT ~ 30 and 
its typical resistance 350 Q. A 250 Q parasitic resistance 
(film residual resistance and CuNi wiring resistance) also 

1 In the literature, another denominator is sometimes used : the 
kinematic viscosity based on normal fluid viscosity and He-II total 
density (0.90. 10 -8 m 2 /s at 1.6K [15]). In our case, it would result 
in Reynolds numbers a decade higher 



contributes to the total electrical resistance. Given such 
conditions (250 < 350), we chose to voltage bias the ther- 
mometer to prevent instabilities which can occur with a 
current bias across positive temperature coefficient ther- 
mistance. The typical voltage bias is a few mV. Facing this 
thermometer, a chromium heating film is deposited on the 
tip of the other beam. The cavity's volume is therefore ge- 
ometrically estimated to be ft ~ h.h.D ~ 1 mm x 1 mm x 
250 fim. Mechanical assembly and electrical connections 
are provided at the other ends of both beams by a stack 
of 3 pieces of silicon wafers. The intermediate one sets the 
gap of the cavity and support on both sides the contact 
pads. The whole geometry is designed to reduce the influ- 
ence on the flow. The micromachining process, assembly 
and packaging of the probe will be described in details 
elsewhere. 

Figure 3 shows the first thermal resonances of a typical 
probe with the fluid at rest. Abscissa is the frequency of 
the applied heating and ordinate is the magnitude of the 
measured temperature oscillation (arbitrary units). The 
n th mode resonance frequency is in good agreement with 
f n = n.c/2D, where c ~ 20 m/s is the second sound veloc- 
ity at 1.6 K [15]. For all the probes tested, the quality fac- 
tor Q of the first resonance is of order 20. It is convenient 
to define £o = n-fi/c.Q as the equivalent bulk dissipation 
by length unit which accounts for the total attenuation on 
the first resonance. 

Figure 3 also shows the resonances with an He flow. 
Compared to the rest situation, the flow causes an extra 
attenuation on the temperature signal. A small fraction 
of it simply results from advection of the thermal wave 
out of the cavity. A simple ballistic model of this effect 
gives the equivalent bulk dissipation ^ a d v = V/2h.c which 
turns out to be one order of magnitude below the signal of 
interest and will be neglected. The remaining of the extra 
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Fig. 3: The first modes of second-sound resonances of the 
probe, with flow (dash line) and without (solid line). 



attenuation is attributed to a bulk dissipation by vortex 
lines ^vld ' this is the signal of interest. Theory predicts 
that a second-order resonance frequency shift is generated 
by vortex dissipation ( [17] and references within) but, in 
our case, it is much smaller than the resonance bandwidth 
and we can neglect this effect. To probe the fluctuations 
of the VLD during experiments, the sensor is operated 
at the first resonance frequency. The measured tempera- 
ture oscillation, called AXb (no flow) and AT (with flow) 
are demodulated by a lock-in amplifier with a 160 fis time 
constant and the oscillation's amplitude is stored for post- 
processing. Basic oscillator theory relates VLD dissipation 
ivLD and the time dependent amplitude AT according to : 



AT 



AT nnh[(Zo+£vLD).D] (1) 

Attenuation ^vld is related to the projected vortex line 
density L± by : 



IVLD 



B.k.L±_ 
4c~~ 



(2) 



where the projected vortex line density L± is defined 



as : 



L 



- I sin 2 O.dl 



(3) 



in which Q is the cavity volume, the summation is per- 
formed along all the vortex line elements dl located inside 
the cavity and is the angle between the line elements 
and the axis of propagation of the second sound waves. 

Heaters in He-II are well known to induce counterflows 
and -for large enough heating- the counterflow can un- 
dergo self-induced turbulent transitions [3-5,18]. Conse- 
quently, a special attention has been dedicated to properly 
choose the heating supplied to the probe, in order to keep 
the probe non-invasive. With the fluid at rest, a tran- 
sition has been evidenced on the temperature oscillation 
for a 30 fiW heating in the Cr film. The corresponding 
heat density, once amplified by the resonance gain Q is 
typical of the turbulence thresholds of counterflows tur- 
bulence [3-5]. During turbulence experiments, we oper- 
ated the probe with driving power ranging up to 10 times 
higher that the above threshold, and down to 5 times be- 
low it. The normalized measured spectrum turns out to 
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Fig. 4: Power spectrum density of the vortex line density L± 
for different mean flow velocities : from bottom to top 0, 0.68, 
0.90 and 1.25 m/s. The straight line is a (-5/3) power law. The 
insert is a f~ 5 ^ 3 compensated spectrum for the 3 different mean 



flows after removal of a 5.10 



Hz white noise floor. 



be independent of the driving heat, indicating that any 
self-induced turbulence transition in this range has a neg- 
ligeable contribution on the signal. It is worth noting that 
Holmes and Van Sciver [12] conducted some second sound 
attenuation measurements with a probe heating density 
more than 4 decades above ours, across a He-II flow at 
similar temperature and velocites. These authors find av- 
erage VLD fully consistent with ours. 

To estimate the signal bandwidth, one needs to consider 
three time constants. First, the geometrically-limited time 
resolution of the probe : the time of flight h/V of or- 
der 1 ms. Second, the time constant of the resonator is 
Qjfx — 0.5 ms, where Q is the quality factor. Third, the 
0.6 ms time constant of the electronic set-up (—3dB cut- 
off for cable/demodulation/acquisition) which has been 
measured without probes, using two frequency mixers to 
mimic the Joule effect frequency doubling and the modula- 
tion of the signal by vortex-line-attenuation. These three 
time scales are comparable and the expected bandwidth 
is therefore DC — 1 kHz. 

Result. — Figure 4 presents the main result of this 
letter : the power spectrum 2 of the vortex line density L± 
for mean flow velocities of 0, 0.68, 0.90 and 1.25 m/s. The 
straight line eyes-guides the —5/3 power law. The insert in 
figure 4 shows the 3 same spectra after subtraction of the 
noise floor, which is fitted from the zero- velocity spectrum 
by a 5.10 15 m~ 4z Hz~ 1 white noise. The low frequency satu- 
ration is consistent with the integral scale plateau of order 
1 cm. In between, over almost one decade, a local power 
exponent of 1.55 ± 0.15 is found. During other runs, this 
scaling has been observed down to the lowest stable mean 



2 The power spectrum density L±(f).L±(f) is normalized such 
that its integral over positive frequencies equals < (Lj_— < Lj_ > 
) 2 >, where < . > denotes time averaging. A Welch algorithm is 
used for spectral averaging. 
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velocity reachable in our set-up (0.3 m/s). The observed 
dependance is compatible with a -5/3 scaling. 

The inverse root of the mean VLD is often seen as a typ- 
ical inter- vortex spacing. This assumes a relative smooth- 
ness of the vortex lines at small scales (negligible "excess- 
length" at scales smaller than inter- vortex spacing). With 
this assumption and assuming isotropy of the tangle, we 
find an inter- vortex spacing of 4 fim for a mean velocity of 
1 m/s. 

We carried several tests to check the nature of the mea- 
sured signal. Among these, we checked that the probe's 
output signal near / wasn't modulated by residual tem- 
perature fluctuations through an hypothetical non-linear 
response of its thermometer. By comparison with the sig- 
nal delivered by a specially designed miniature velocity 
probe [16], we also checked that the second sound probe 
wasn't simply measuring velocity, for exemple through the 
vorticity generated in the thin boundary layers. We also 
checked that the measured signal wasn't history depen- 
dent, as it sometimes happens with some superfluid sen- 
sors due to trapped vortex lines [19]. 

Interpretation and Conclusion. — The interpreta- 
tion of the observed spectrum is still open. As a con- 
clusion, we discuss two approaches of it. Due to a lack 
of theoretical prediction and numerical observation of the 
projected VLD L± power spectrum, we tried in both case 
to relate L±_ to a classical counterpart. 

A first approach would amount to relate the VLD to 
the kinetic energy, or more precisely L± to second order 
statistics of the velocity components. Indeed, fourth order 
velocity statistics, such as power spectrum of the second 
order velocity components, are known to produce —5/3 
power laws ( [21-23]). In a fully random tangle, the VLD 
is likely to be proportional to the kinetic energy ( [4] page 
50-51), at least up to a log correction and assuming that 
the vortices are straight enough at scales smaller than the 
inter-vortex spacing. This linear relation can account for 
a —5/3 spectrum. Unfortunately, as pointed to us by J. 
Vinen, this simple relation probably no longer holds in a 
Kolmogorov-like tangle : the partial polarization of vor- 
tices should be strong enough to invalidate the random 
tangle hypothesis. A more detailed modeling is therefore 
necessary here. 

Alternatively, the VLD can also be seen as the quantum 
analog of the enstrophy (vorticity square) or the absolute 
value of the vorticity in classical fluids [20] . Enstrophy in 
classical turbulence follows a +1/3 power law spectrum. 
Unfortunately we couldn't find any reports on the corre- 
sponding spectrum for the projected enstrophy (projected- 
vorticity square) nor for the absolute value of the pro- 
jected vorticity. Nevertheless, with the sweeping theory 
in mind [21,22], we can conjecture that the inertial-range 
dynamics of the polarization of the vortices is controled 
by energy containing eddies, which would result in the ob- 
served —5/3 power law. This picture has some analogies 
with passive scalar fluctuations, which also have a similar 



scaling. 

* * * 

The idea of an in-flow probe was proposed by P. Tabel- 
ing and macroscopic protoypes were done with the help 
of J. Maurer at the Ecole Normale Superieure. Present 
work results from a collaboration between three laborato- 
ries. The ESPCI designed the emitter and receiver parts 
of the probes, and made the Al thermometry. The ESIEE 
micro-machined the emitter and receiver. The CRTBT 
completed the probe design / machining and made the 
experiment and analysis. 
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